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■ A new class of automata networks is denned. Their evolution rules are determined by 

a probability measure p on the set of all integers Z and an indicator function I a on 
the interval [0, 1]. It is shown that any cellular automaton rule can be represented by 
a (nonunique) rule formulated in terms of a pair (p, I a)- This new class of automata 
networks contains discrete systems which are not cellular automata. For a given p, 
a metric can be defined on the space of all rules which induces a metric on the space 
' of all cellular automata rules. 
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Systems which consist of a large number of simple identical elements evolving 
in time according to simple rules often exhibit a complex behavior as a result 
of the cooperative effect of their components. Automata network are models 
of such systems. They consist [9] of a graph with a discrete variable at each 
vertex which evolves in discrete time steps according to a definite rule in- 
volving the values of neighbouring vertex variables. The vertex variables may 
be updated sequentially or synchronously. More formally, automata networks 
may be defined as follows: 

Let G = (V, E) be a graph, where V is a set of vertices and E a set of edges. 
Each edge joins two vertices not necessarily distinct. An automata network, 
defined on V, is a triple (G, Q, {fi\i G V}), where G is a graph on V, Q a finite 
set of states and fi'-Q^ — > Q a mapping, called the transition rule associated 
to vertex i. Ui = {j G V\{j, i} G E} is the neighbourhood of i, i.e., the set of 
vertices connected to i, and \Ui\ denotes the number of vertices belonging to 
Ui. If \Ui\ is finite, the rule is local. For the automata networks considered in 
this paper, \Ui\ is infinite, however it will be shown that, if certain conditions 
are fulfilled, even if the number of vertices in the neighborhood of a vertex is 
infinite, the rule may be local. Cellular automata (CA) are particular automata 
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networks in which the set of vertices V is the set of all integers Z , where d is 
called the space dimensionality, and the rule is translationally invariant, local 
and applied synchronously. 

Since the class of automata networks to be described is a generalization of 
CA, to fix the notations, we define CA as follows: Let s : Z x N i— > {0, 1} be 
a function that satisfies the equation 

(Vi G Z) (Vt G N) s(i, t + l) = f(s(i -r,t), 

s(i-r + l,t),...,s(i + r,t)), (1) 

where N is the set of nonnegative integers and Z the set of all integers. Such 
a discrete dynamical system is a two-state one-dimensional CA. The mapping 
/ : {0, l} 2r+1 — > {0,1} is the rule, and the positive integer r is the radius 
of the rule. More general rules could be site- and/or time-dependent and/or 
involve different left and right radii. The function St : % i— > s(i, t) is the state of 
the CA at time t. S = {0, 1} Z is the state space. An element of the state space 
is also called a configuration. Since the state S t +i at time t + 1 is entirely 
determined by the state St at time t and the rule /, there exists a unique 
mapping Ff : S — > S such that St+i = Ff(St). Ff, which is the evolution 
operator, is also referred to as the CA rule. 

CA have been widely used to model complex systems in which the local char- 
acter of the rule plays an essential role [8,16,13,10,3]. Nevertheless, not all 
complex systems exhibit purely local interactions. Consider, for example, dif- 
fusion of innovations in a social system. For a given individual, deciding to 
buy, say, her first computer is based on the word of mouth from, more or 
less, closed neighbors as her family members, friends and colleagues as well as 
mass-media communications. It is clear that a CA rule cannot correctly model 
this type of system. The new class of automata networks which described in 
this paper might be useful to model such systems. 

The evolution operator F p a of our new class of automata networks is defined 
in terms of a probability measure p on Z and an indicator function I a on [0, 1]. 
Then, we shall show that, for any CA rule / , we can find such an evolution 
operator, that is, a measure p and an indicator I a, that emulates the CA rule. 

Let S t be the state of the system at time t, that is, St : i i— > s(i, t), and put 

oo 

a(i,t)= ]T s(i + n,t)p(n), (2) 

n=— oo 

where p is a given probability measure on Z, that is, a nonnegative function 
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on the set of all integers such that 

oo 

E PW = L (3) 



For all i 6 Z and t G N, cr(i,t) G [0, 1]. The state S t+1 of the system at time 
t + 1 is then determined by the function 

i^s(i,t + l)=I A (a(i,t))=I A ( £ s(i + n,t)p(n)J , (4) 



where Ja is a given indicator function on [0, 1], that is, a function such that, 
for all x G [0, 1], 

'l,ifxGAc[0,l], 
Ia(x) = { (5) 
0, otherwise. 



Since the state S t+1 at time t+ 1 is entirely determined by the state S t at time 
£ and the measure p and the subset A of [0, 1], there exists a unique mapping 
F PtA : 5 — > S such that 

«St+i = F P ,A{St)- 
F p< a will be referred to as the evolution operator. 

For a given probability measure p and a given indicator function I a, the op- 
erator F Pi a is entirely determined. That is, the word "probability measure" 
should not be misleading: the evolution operator F p ^a is deterministic. Prob- 
abilistic evolution operators could also be defined. If, for instance, we replace 
Ia by XIa, where X is a Bernoulli random variable taking values on {0, 1}, 
the resulting evolution operator would be probabilistic. 

Theorem 1 Let Ff be the evolution operator on S associated to a CA rule f; 
then, there exists an evolution operator F Pt A such that, for any configuration 
x G S, F P)A {x) = F f (x). 

To simplify the proof, we only consider translation-invariant symmetric CA 
rules, that is, rules / such that, for all i e Z and all (2r + l)-block B(i,r) = 
{x(i — r),x(i — r + 1), . . . , x(i + r)} (r > 0), 

f(x(i — r),x(i — r + 1), . . . , x(i + r)) = 

f(x(i + r),x(i + r — 1), . . . , x(i — r)). (6) 
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The generalization to nonsymmetric rules is straightforward. Since / is sym- 
metric, we shall verify that we may assume that p is even, that is, for all 
n G Z, p(—n) = p{n). The set of configurations with prescribed values at a 
finite number of sites is called a cylinder set. To each (2r + l)-block B(i, r) cor- 
responds a cylinder set denoted C(i,r). Since the rule is translation- invariant, 
it is sufficient to consider i — 0. 

For any configuration i:jih x(n) belonging to the cylinder set C(0, r), the 
set of all numbers £(C(0, r)) defined by 

oo 

e(C(0,r))= £ s(n)p(n) (7) 

n=— oo 

belongs to the subinterval [£ m in(C(0, r)), £ max ((7(0, r))] of [0,1] (called a C- 
interval in what follows) such that 

r r 

6nin(C(0,r))= £ x(ra)p(ra) = x(0)p(0) + £ (x(-ra) + x(n))p(n) (8) 

n=— r n=l 

and 

oo 

£max(C(0,r))=£ m in(C(0,r)) + 2 £ p ( n)> (9) 

n=r+l 

where we have taken into account that p(—n) = p(n). 

Since, for n ^ 0, £ m in(C(0, r)) depends on x(—n) + x(n) and not on x(—n) 
and x(n) separately, there are only 2 x 3 r different ^ m in(C(0, r)) (and as many 
C-intervals) . It is, therefore, more convenient to label them £ m i n (z/), where v 
is an integer in {0, 1, . . . , 2 x 3 r — 1} defined by 

r 

u = x(0)3 r + J2«-k)+x(k))3 r ~ k . (10) 
k=i 

The corresponding £ m i n (v) is then given by 

= x(0)p(0) + £ (x(-k) + x(k))p(k). (11) 
k=i 

If we want to define 2 2xy different operators F Pi a representing the 2 2x3r dif- 
ferent symmetric CA rules / of radius r, we have first to find a probability 
measure p such that the 2 x 3 r C-intervals are disjoint. Then, according to the 
rule / to be represented, we will choose a subset A of [0, 1] such that some of 
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these C-intervals are strictly included in A whereas the others have an empty 
intersection with A. The only problem is, therefore, to find the conditions to 
be satisfied by p. 

The sequence {{ m in(^) | v = 0, 1, . . . , 2 x 3 r — 1} is totally ordered for increasing 
values of v if, and only if, the r + 1 conditions 

p(r -k)> 2p(r - k + 1) + 2p(r - k + 2) • • • + 2p(r) (12) 

(fc = 0,l,...,r) 

are satisfied, with the convention: p{m) = if m < 0. 

It is easy to verify that, for v — 0, 1, . . . , 2 x 3 r — 1, the difference £ m i n (^ + 
1) — £ m i n (^) is always of the form 

p(r — k) — 2p(r - k + 1) - 2p(r - k + 2) 2p(r) (13) 

(fc = 0,l,...,r). 

Therefore, if all these differences are greater than 

oo 

W(C(0,r))-Uin(C(0,r)) = 2 £ 

n=r+l 

all the intervals 

oo 

^min(^)) £min 

(z/) + 2 ^ 

n=r+l 

will be disjoint. For this to be the case, the r + 1 conditions 

oo 

p(r-k)>2 P( n ) (fc = 0,l,...,r). (14) 

n=r— fc+1 

should be satisfied. Since p satisfies (3), Conditions (14) are identical to Con- 
ditions (12). 



2 Examples 

Here are two examples of probability measures p. 
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Example 1. If, for all n e Z, 

p(n) = tanh(|A) exp(— A|n|), (15) 

where A > 0, it is easy to verify that Conditions (12) are satisfied for all 
positive values of r if A > log 3. That is, choosing a subset A of [0, 1] we can 
represent, in this case, any symmetric CA rule. For a given CA rule p and A 
are clearly not unique. 

If we take A = 1.2, the corresponding C-intervals are approximately 

[£mm(0),£max(0)] = [0,0.1395] 
[&nin(lUmax(l)] = [0.1618, 0.3012] 

[£min(2),£nax(2)] = [0.3235, 0.4630] 
[£min(3),£nax(3)] = [0.5370, 0.6765] 
[6nin(4), £max(4)] = [0.6988, 0.8382] 

[6nin(5),£nax(5)] = [0.8606,1]. 

And to obtain Rule 18 [15] defined by 



f(x 1 ,x 2 ,x 3 ) = 



1, if (x 1 ,x 2 ,x 3 ) = (0,0,1) or (1,0,0), 
0, otherwise, 



we may choose A = [0.15, 0.31]. 
Example 2. If, for all n G Z, 

where a > 1 and ( is the zeta function, then, for a radius equal to ro, there 
exist a threshold value ao such that Conditions (12) are satisfied for all r < r 
if a > a where a is the solution of the equation 



(l + r )" ao -2 J2 (l + n)- ao = 0. 

n=ro+l 

The table below give some approximate threshold values. 



r 1 2 3 10 20 100 

a 2.18529 3.28994 4.39061 5.49023 13.1824 24.169 112.058 
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a a is approximately a linear function of tq. 



For a = 3.5, the corresponding C-intervals are approximately 



[£min(0Umax(0)] = [0,0612] 
[£min(lUmax(l)] = [0.0705, 0.1317] 
[£min(2),£nax(2)] = [0.1410, 0.2022] 
[£min(3Umax(3)] = [0.7978, 0.8590] 

[£ min (4),£ max (4)] = [0.8683, 0.9295] 

[Un(5),emax(5)] = [0.9388,1]. 

And to obtain Rule 18, we may choose, in this case, A = [0.062,0.133]. 

If some C-intervals are not disjoint but are either strictly included in A or 
have an empty intersection with A, certain CA rules could be represented, 
but not all CA rules of a given radius. To represent nonsymmetric CA rules 
we should consider noneven functions p. 

An operator F PjA will not, in general, represent a CA rule in the following two 

cases. 

(i) If p is such that for a given r the C-intervals are disjoint but A has a 
nonempty intersection with some C intervals without, however, strictly 
including these intervals. 

(ii) If p is such that for a given r the C-intervals are not disjoint. 

In the first case, the spatio-temporal pattern will look like the pattern of a CA 
whose evolution is governed by a mixed rule. But, a mixed rule being local, 
the velocity of propagation of a local perturbation cannot be greater than r, 
whereas in this case the value of the a-variable defined by (2) depends on the 
whole configuration and not only on a finite block. Note that the nonlocal 
character of the evolution operator F Pj a comes essentially from the definition 
of the indicator function since if we consider a measure p whose moments are 
all defined (as in Example 1) the range of the interaction between the different 
sites is clearly finite whatever the precise meaning we give to this range. As 
a consequence of the nonlocal character of the rule in this case the velocity 
of propagation of a local perturbation may be time-dependent and grow to 
infinity. 

In the second case, the spatio-temporal pattern may look like a CA pattern, 
but no finite radius can be defined. As an example, consider the exponential 
probability measure, and choose A = 1. The C-intervals are not disjoint and 
the resulting evolution operator is not a CA rule. For two different initial 
configurations, the corresponding spatio-temporal patterns of such an operator 
are represented in Figures la and lb. While Pattern la is reminiscent of a 
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Fig. 1. Spatiotemporal pattern for A = [0.197,0.600] and A = 1 for two different 
initial configurations. 



class-2 CA, Pattern lb looks more like class-3. Wolfram [15] classification is, 
therefore, not relevant for rules of this type (failure of Wolfram classification 
has already been noticed even for standard CA, as reported in [11] and [7]). 

CA have been widely used to model systems, say in epidemiology [1,2,4,5] or 
in ecology [5], in which the local character of the interaction is an essential 
feature. But, there are cases in which the interactions are both local and 
nonlocal as for the diffusion of innovations [12,14] which can be influenced by 
two types of communication channels: mass media and interpersonal. The new 
class of discrete dynamical systems we have defined may be well adapted to 
model such a situation. If, for instance, we consider a probability measure p 
such that p{n) is proportional to 



where Ai, A 2 and a are positive constants, and n is a positive integer. Ai and 
A2 are inversely proportional to the range of the nearby and distant sources 
of information, a measures the relative weight of the distant source compared 
to the nearby one and n is the location of the distant source. A model of 
diffusion of innovations employing this idea is currently investigated and a 
detailed report will be published elsewhere. 
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3 Metric 



For a given probability measure p, let F be the probability distribution of the 
random variable E : S — > [0, 1] defined by 

oo 

E = ]T X(n)p(n), (18) 

n=— oo 

where {A(n) n G Z} is a doubly infinite sequence of equally distributed 
Bernoulli random variables such that 

(Vn G Z) P(X(n) = 0) = P{X{n) = 1) = ±. (19) 

Theorem 2 If \ > log(3), i/ie random variable E is singular. 

That is, for A > log 3, the distribution function F of E is continuous but 
has no density. The derivative of F is almost everywhere equal to zero, and 
F increases on a Cantor-like set of zero Lebesgue measure (see Figures 2a 
and 2b). This implies that, if an evolution operator is selected at random 
by selecting a subset A of [0, 1] at random, then, with probability one, the 
corresponding rule will be a CA rule. 

As for the usual triadic Cantor set, the set on which E is defined can be 
constructed by removing step by step "forbidden intervals." 

1st step 

oo 

<7^]2 5>(n),p(0) 

n=l 

2nd step 

^]2&W,P(1)' 

n=2 

oo 

U]p(l) + 2$>(n),2p(l)' 

n=2 

oo 

U]p(0) + 2]Tp(n),p(0)+p(l)' 

n=2 

oo 

U]p(0) + p(l) + 2 £ p(n),p(0) + 2p(l) 

n=2 

etc. For a given radius r, the sum m(r) of the Lebesgue measures of the 
removed intervals is 
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a) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



a 



b) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



a 



Fig. 2. Distribution function for a) A = 1.3 and b) A = 2 . 

2x3 r -l 

m ( r ) = (£min(» - £max(^ ~ 1)) 

u=l 

= (p(0) + 2p(l) + ■■■ + 2p(r)) 

CO 

+ (2x3 r -l)2 Yl P(k)- 

k=r+l 
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Since 

oo 

hm (2 x 3 r - 1) 2 ]T p(fc)=0. (21) 

fc=r+l 



it follows that 

lim m(r) = 1. (22) 



This result could be made more intuitive in the following way. If we put 
A = log b, 

E ^ m+ g^2My (23) 



Since, for all positive integers n, the sum of the Bernoulli random variables 
X(—n) and X(n) is either equal to 0, 1 or 2, the summation over n in (23) 
represents a random variable whose values are expressed in base b {b is not 
necessarily an integer). If b > 3, the set of all possible values of £ is, therefore, 
a Cantor-like set. 

For 6 = 3, the random variable £ is given by 



X(—n) + X{n) being a binomial random variable such that 



P(X(-n)+X(n) = 0) = i 
P(X(-n)+X(n) = l) = \ 
P(X(-n)+X(n) = 2) = \ 

£ is absolutely continuous but not uniformly distributed on [0, 1] (see Figure 
3), and the set of all possible values of £ has a nontrivial Hausdorff dimension 
dn, which, in this very particular case, can be exactly determined (see, for 
example, [6]). It is found that 

d H = --i- = 0.946395.... (25) 
21og3 V ; 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



a 

Fig. 3. Distribution function for A = log 3. 

For a fixed probability measure p we define the distance between two evolution 
operators F PjAl and F p ^ 2 by 

l 

d(F PtAl ,F PtAl ) = J \I Al (a)-I A2 (a)\dF(a) = J dF(a), (26) 

A 1 AA 2 

where A1AA2 denotes the symmetrical difference between A\ and A2, that is, 
the set of all numbers a G [0, 1] that belong to A± XOR A 2 . 

Note that d(F P:Al ,F PjAl ) = does not imply A\ = A 2 . This means that two 
evolution operators may be different but the corresponding dynamics will be 
identical. Therefore, we could, instead of considering evolution operators we 
could consider classes of evolution operators. Two evolution operators would 
belong to the same class if their distance is equal to zero. This notion of 
distance may be useful to study sequences of, say, cellular automaton rules 
and study if such sequences converge in the topology defined by this metric. 
This topic will be developed elsewhere. 

Theorem 3 If X > log(3), the distance between two CA rules does not depend 
upon X. 

This result follows from the fact that the subset A, characterizing a CA rule, 
is a reunion of disjoint intervals whose end points belong to removed intervals. 
It can even be proved that the distance will not depend upon the form of p 
provided that p is such that Conditions (12) are satisfied for all values of r less 
or equal than the larger radius of the rules under consideration. For example, 
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the distance between Rule 18 (denned above) and Rule 22 denned by 



f(x!,x 2 ,x 3 ) 



1, if (x 1 ,x 2 , x 3 ) = (0, 0, 1) or (1, 0, 0) or (0, 1,0), 
0, otherwise, 



is 



d(F 18 ,F 22 ) = J dF(a) = ^ 



(27) 



{0,1,0} 



It is equal to the probability of selecting {0, 1, 0} among the 8 equally probable 
blocks of radius one. 



4 Conclusion 

We have described a new class of automata networks whose evolution opera- 
tors are defined in terms of a probability measure p and an indicator function 
I a, where A is a subset of the interval [0,1]. This class contains all cellular 
automaton rules, and other rules which have no finite radius. For a given p, 
we have defined a metric on the space of these new rules which induces a 
metric on the space of all CA rules that does not depend upon the particular 
p. In the case of evolution operators which are not CA rules, the Wolfram 
classification into 4 classes does not seem to be relevant since we may ob- 
tain spatio-temporal patterns reminiscent of two different classes by changing 
the initial configuration. The parameters characterizing the evolution opera- 
tors of these new automata networks can be modified continuously, and, as 
a consequence, starting from a given initial configuration, we can, say, pass 
continuously from a class-3 to a class-2 CA. Finally these new rules could 
be useful to model systems in which the interactions between the different 
elements are short- and long-range, like in the diffusion of innovations. 
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